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Abstract 



The theoretical basis for the prediction of anisotropies in the cosmic microwave background 
, is very well developed. Very low amplitude density and temperature perturbations produce 

lO ' small gravitational effects, leading to an anisotropy that is a combination of temperature 

fluctuations at the surface of last scattering and gravitational redshifts both at last scattering 
and along the path to the observer. All of the primary anisotropy can be handled by linear 
perturbation theory, which allows a very accurate calculation of the predicted anisotropy 
from different models of the Universe. 
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Pi 1.1 Introduction 

' The first predictions of the anisotropy of the cosmic microwave background (CMB) 

were published shortly after the CMB was discovered by Penzias & Wilson (1965). Sachs & 
Wolfe (1967) calculated the anisotropies due to gravitational potential fluctuations produced 
by density perturbations (Figure [T~T1 . Because the density perturbations are given by the 
second derivative of the gravitational potential fluctuation in Poisson's equation, the Sachs- 
Wolfe effect dominates the temperature fluctuations at large scales or low spherical harmonic 
index I, Sachs & Wolfe predicted AT /T « 10~ 2 at large scales. This prediction, which failed 
by a factor of 10 3 , is based on correct physics but incorrect input assumptions: prior to the 
discovery of the CMB no one knew how uniform the Universe was on large scales. 

Silk (1968) computed the density perturbations needed at the recombination epoch at 
Z Ri 10 3 in order to produce galaxies, and predicted AT /T «3x 10~ 4 on arcminute scales. 
Silk (1967) calculated the damping of waves that were partially optically thick during re- 
combination. This process, known as "Silk damping," greatly reduces the CMB anisotropy 
for small angular scales. 

Observations by Conklin (1969) and then Henry (1971) showed that there was a dipole 
anisotropy in the CMB corresponding to the motion of the Solar System with respect to the 
average velocity of the observable Universe. There is a discussion of the dipole observations 
and their interpretation in Peebles (1971) that is still valid today, except that what was then 
a "tentative" dipole is now known to better than 1% accuracy, after a string of improved 
measurements starting with Corey & Wilkinson (1976) and ending with the COBE DMR 
(Bennett et al. 1996). 

Peebles & Yu (1970) calculated the baryonic oscillations resulting from interactions be- 
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tween photons and hydrogen in the early Universe, and also independently introduced the 
Harrison-Zel'dovich spectrum (Harrison 1970; Zel'dovich 1972). Pebbles & Yu predicted 
AT/T w 1.5 x 10" 4 on 1' scales and AT/T w 1.7 x 10" 3 on 7' scales. 

Wilson & Silk (1981) further developed the theory of photon and matter interaction by 
scattering and gravity, and predicted AT/T = 100 ,uK for a single subtracted experiment 
with a 7° throw and with a 7° beam like COBE. Of course, when COBE was launched in 
1989 it actually observed a much smaller anisotropy. 

These early predictions of a large anisotropy were greatly modified by the addition of dark 
matter to the recipe for the cosmos. Observational upper limits on small-scale anisotropies 
had reached AT /T < 4 x 10~ 5 on 1 .'5 scales (Uson & Wilkinson 1982), which was consider- 
ably less than the predictions from universes with just baryons and photons. Peebles (1982) 
computed the anisotropy expected in a universe "dominated by massive, weakly interacting 
particles" — in other words cold dark matter (CDM), although this paper predated the use 
of "cold dark matter." 

Bond & Efstathiou (1987) calculated the correlation function of the CMB anisotropy, 
C{6), and also the angular power spectrum, Q, in the CDM cosmology. This paper contains 
one of the first plots showing £(£+ 1)Q vs. £, with peaks originally called the "Doppler" 
peaks but more properly called "acoustic peaks." This paper solved the Boltzmann equa- 
tion describing the evolution of the photon distribution functions. Several authors developed 
these "Boltzmann codes," but the calculation of the angular power spectrum up to high I 
was very slow. These codes described the conversion of inhomogeneity at the last-scattering 
surface into anisotropy on the observed sky by a set of differential equations evolving the co- 
efficients of a Legendre polynomial expansion of the radiation intensity. Since the Universe 
is almost completely transparent after recombination, a ray-tracing approach is much more 
efficient. This great step in efficiency was implemented in the CMB FAST code by Seljak & 
Zaldarriaga(1996). 

Hu & Dodelson (2002) give a recent review of CMB anisotropies, which includes a very 
good tutorial on the theory of AT/T. 



2 



E. L. Wright 



1.2 Results 

A simple analysis of cosmological perturbations can be obtained from a consider- 
ation of the Newtonian approximation to a homogeneous and isotropic universe. Consider 
a test particle at radius R from an arbitrary center. Because the model is homogeneous the 
choice of center does not matter. The evolution of the velocity of the test particle is given by 
the energy equation 

v 2 GM 

-J^ m+ —. (1.1) 

If the total energy E tot is positive, the Universe will expand forever since M, the mass (plus 
energy) enclosed within R, is positive, G is positive, and R is positive. In the absence of a 
cosmological constant or "dark energy," the expansion of the Universe will stop, leading to 
a recollapse if E, , is negative. But this simple connection between E tot and the fate of the 
Universe is broken in the presence of a vacuum energy density. The mass M is proportional 
to R 3 because the Universe is homogeneous and the Hubble velocity v is given by v = HR. 
Thus E,„, oc R 2 . 

We can find the total energy by plugging in the velocity vq = HqRq and the density po in 
the Universe now. This gives 

_ (HqRqc) 2 4nGp R 2 _ (H Q R U ) 2 ( p Q \ 
"" 2 3 2 \ Pcr J ' 

with the critical density at time to being p„„ = 3//q/(87tG). We define the ratio of density to 
critical density as £1 = p/p cr it- This includes all forms of matter and energy. fl m will be 
used to refer to the matter density. 

From Equation ll.il we can compute the time variation of il. Let 

, 2GM , , SirGpR 2 

2E tot = v 2 -——=H 2 R 2 £—= const. (1.3) 

R j 

If we divide this equation by %-irGpR 2 /3 we get 

m " 1 = ^ = ^-1. (1.4) 



SirGp pR 2 

Thus fT 1 — 1 oc (pR 2 )~ l . When p declines with expansion at a rate faster than R~ 2 then the 
deviation of from unity grows with expansion. This is the situation during the matter- 
dominated epoch with p oc R~ 3 , so il" 1 - 1 oc R. During the radiation-dominated epoch p oc 
R~ A , so JT 1 - 1 oc R 2 . For f2 to be within 0.9 and 1.1, O needed to be between 0.999 and 
1.001 at the epoch of recombination, and within 10~ 15 of unity during nucleosynthesis. This 
fine-tuning problem is an aspect of the "flatness-oldness" problem in cosmology. 

Inflation produces such a huge expansion that quantum fluctuations on the microscopic 
scale can grow to be larger than the observable Universe. These perturbations can be the 
seeds of structure formation and also will create the anisotropies seen by COBE for spherical 
harmonic indices i > 2. For perturbations that are larger than ~ c s t (or ~ c s /H) we can 
ignore pressure gradients, since pressure gradients produce sound waves that are not able to 
cross the perturbation in a Hubble time. In the absence of pressure gradients, the density 
perturbation will evolve in the same way that a homogeneous universe does, and we can use 
the equation 
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pa 2 1^ = const, (1.5) 

the assumption that Q m 1 for early times, and Ap <C p as indicated by the smallness of the 
Ar's seen by COBE, to derive 

-pa 2 ( i-1 J ~ pcrna 2 Afl w Apa 2 = const. (1.6) 



Hence, 



GAM 4tt GAppjaL) 3 1 A Po 2 

A</>= — — = - = - (H L) , (1.7) 

« 3 aL 2 p crir 

where L is the comoving size of the perturbation. This is independent of the scale factor, so 
it does not change due to the expansion of the Universe. 

During inflation (Guth 2003), the Universe is approximately in a steady state with con- 
stant H. Thus, the magnitude of A<j) for perturbations with physical scale c/H will be the 
same for all times during the inflationary epoch. But since this constant physical scale is 
aL, and the scale factor a changes by more than 30 orders of magnitude during inflation, 
this means that the magnitude of Acj) will be the same over 30 decades of comoving scale 
L. Thus, we get a strong prediction that A<j> will be the same on all observable scales from 
c/Hq down to the scale that is no longer always larger than the sound speed horizon. This 
means that 

— ocZT 2 , (1.8) 
P 

so the Universe becomes extremely homogeneous on large scales even though it is quite 
inhomogeneous on small scales. 

This behavior of A<p being independent of scale is called equal power on all scales. It was 
originally predicted by Harrison (1970), Zel'dovich (1972), and Peebles & Yu (1970) based 
on a very simple argument: there is no scale length provided by the early Universe, and thus 
the perturbations should be scale-free — a power law. Therefore A(f> oc LT. The gravitational 
potential divided by c 2 is a component of the metric, and if it gets comparable to unity then 
wild things happen. If m < then A<p gets large for small L, and many black holes would 
form. But we observe that this did not happen. Therefore m > 0. But if m > then A</> gets 
large on large scales, and the Universe would be grossly inhomogeneous. But we observe 
that this is not the case, so m < 0. Combining both results requires that m = 0, which is a 
scale-invariant perturbation power spectrum. This particular power-law power spectrum is 
called the Harrison-Zel'dovich spectrum. It was expected that the primordial perturbations 
should follow a Harrison-Zel'dovich spectrum because all other answers were wrong, but 
the inflationary scenario provides a good mechanism for producing a Harrison-Zel'dovich 
spectrum. 

Sachs & Wolfe (1967) show that a gravitational potential perturbation produces an aniso- 
tropy of the CMB with magnitude 

Ar 1 A6 

- = 3^ d.9) 

where A<f> is evaluated at the intersection of the line-of-sight and the surface of last scattering 
(or recombination at z ~ 1 100). The (1 /3) factor arises because clocks run faster by a factor 
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Fig. 1.2. Left: A plane wave on the last-scattering hypersurface. Right: The spherical 
intersection with our past light cone is shown. 



(1 +<fi/c 2 ) in a gravitational potential, and we can consider the expansion of the Universe 
to be a clock. Since the scale factor is varying as a oc f 2 ' 3 at recombination, the faster 
expansion leads to a decreased temperature by AT /T = -(2/3)A0/c 2 , which, when added 
to the normal gravitational redshift AT /T = A(p/c 2 yields the (1 /3) factor above. This is 
an illustration of the "gauge" problem in calculating perturbations in general relativity. The 
expected variation of the density contrast as the square of the scale factor for scales larger 
than the horizon in the radiation-dominated epoch is only obtained after allowance is made 
for the effect of the potential on the time. For a plane wave with wavenumber k we have 
-k 2 A(j) = 4irGAp, or 



^=-W) 2 ^, (1.10) 
2 p crit 



so when p w p crit at recombination, the Sachs-Wolfe effect exceeds the physical temperature 
fluctuation AT /T = (l/4)Ap/p by a factor of 2(H/ck) 2 if fluctuations are adiabatic (all 
component number densities varying by the same factor). 

In addition to the physical temperature fluctuation and the gravitational potential fluc- 
tuation, there is a Doppler shift term. When the baryon fluid has a density contrast given 
by 

5b(x,t)= = 5i,exp[ik(x—c s t)], (111) 
Pb 

where c s is the sound speed, then 

— — =-ikc s S h p h = -p b V -v = ikvp h . (1-12) 
at 

As a result the velocity perturbation is given by v = -c s Si,. But the sound speed is given 
by c s = y/dP/dp = ^(4/3)p 7 c 2 /(3pb + 4py) w c/\f?>, since p 7 > p b at recombination (z = 
1100). But the photon density is only slightly higher than the baryon density at recom- 
bination so the sound speed is about 20% smaller than c/y/3. The Doppler shift term in 
the anisotropy is given by AT /T = v/c, as expected. This results in AT /T slightly less 
than Sb/V^, which is nearly \/3 larger than the physical temperature fluctuation given by 
AT/T = 6 b /3. 
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Single k=50 plane wave k=50 velocity term 




Fig. 1.3. Left: The scalar density and potential perturbation. Right: The vector velocity 
perturbation. 



These plane-wave calculations need to be projected onto the sphere that is the intersection 
of our past light cone and the hypersurface corresponding to the time of recombination. Fig- 
ure ll.2l shows a plane wave on these surfaces. The scalar density and potential perturbations 
produce a different pattern on the observed sky than the vector velocity perturbation. Figure 
ll.3l shows these patterns on the sky for a plane wave with IcRls = 50, where Rls is the radius 
of the last-scattering surface. The contribution of the velocity term is multiplied by cos 9, 
and since the RMS of this over the sphere is yj 1 /3, the RMS contribution of the velocity 
term almost equals the RMS contribution from the density term since the speed of sound is 
almost c/V3. 

The anisotropy is usually expanded in spherical harmonics: 

^Y^- = ^2^2a£mYem(h). (1.13) 

e m=-e 

Because the Universe is approximately isotropic the probability densities for all the different 
m's at a given I are identical. Furthermore, the expected value of AT(h) is obviously zero, 
and thus the expected values of the ag m 's is zero. But the variance of the a&„'s is a measurable 
function of I, defined as 

d = (\a em \ 2 ). (1.14) 

Note that in this normalization Q and ag m are dimensionless. The harmonic index £ as- 
sociated with an angular scale 9 is given by I i=s 180°/$, but the total number of spherical 
harmonics contributing to the anisotropy power at angular scale 9 is given by A£ « I times 
2£+l. Thus to have equal power on all scales one needs to have approximately Q oc £~ 2 . 
Given that the square of the angular momentum operator is actually £(£ + 1), it is not sur- 
prising that the actual angular power spectrum of the CMB predicted by "equal power on all 
scales" is 

r _ 4 ^(6 2 ) 6 r11 <, 
Q "^T£(ZTT)' ( } 

where (Q 2 ) or Q 2 m ^ PS is the expected variance of the I = 2 component of the sky, which must 
be divided by Tq because the ai, n 's are defined to be dimensionless. The "47r" term arises 
because the mean of \Y^ n \ 2 is l/(47r), so the | a £ m | 2 ' s must be 4ir times larger to compensate. 
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Fig. 1.4. Left: Mean £ is plotted vs. the wavenumber kRis- The light solid line shows £ = 
kRis, while the solid line shows £ = (w/4)kRis, and the dashed line shows £ = (in / \6)kRu. 
Right: The angular power spectra for single k skies is plotted for kRis = 25, 50, and 100 
(from bottom to top). 

Finally, the quadrupole has 5 components, while Q is the variance of a single component, 
giving the "5" in the denominator. The COBE DMR experiment determined y/ (Q 2 ) = 18 /iK, 
and that the Q's from £ = 2 to £ = 20 were consistent with Eauation ll.151 
The other common way of describing the anisotropy is in terms of 



Note these definitions give AT 2 2 = 2.4(Q 2 ). Therefore, the COBE normalized Harrison- 
Zel'dovich spectrum has AT e 2 = 2.4 x 18 2 = 778 piK 2 for £ < 20. 

It is important to realize that the relationship between the wavenumber k and the spherical 
harmonic index £ is not a simple £ = kRis- Figure [T31 shows that while £ = kR^ at the 
"equator" the poles have lower I. In fact, if /j, = cos 6, where 9 is the angle between the wave 
vector and the line-of-sight, then the "local t' is given by kRis-^l — fi 2 . The average of this 
over the sphere is (£) = (ir/4)kRLs- For the velocity term the power goes to zero when /i = 
on the equator, so the average £ is smaller, (£) = {3ir/l6)kR LS , and the distribution of power 
over £ lacks the sharp cusp at £ = kRi$. As a result the velocity term, while contributing about 
60% as much to the RMS anisotropy as the density term, does not contribute this much to 
the peak structure in the angular power spectrum. Thus the old nomenclature of "Doppler" 
peaks was not appropriate, and the new usage of "acoustic" peaks is more correct. Figure 
1 1.41 shows the angular power spectrum from single k skies for both the density and velocity 
terms for several values of k, and a graph of the variance-weighted mean £ vs. kRis- These 
curves were computed numerically but have the expected forms given by the spherical Bessel 
function jg for the density term and f e for the velocity term. 

Seljak (1994) considered a simple model in which the photons and baryons are locked 
together before recombination, and completely noninteracting after recombination. Thus the 
opacity went from infinity to zero instantaneously. Prior to recombination there were two 
fluids, the photon-baryon fluid and the CDM fluid, which interacted only gravitationally. 
The baryon-photon fluid has a sound speed of about cj \/3 while the dark matter fluid has 
a sound speed of zero. Figure 1131 shows a conformal spacetime diagrams with a traveling 
wave in the baryon-photon fluid and the stationary wave in the CDM. The CDM dominates 




'£{£+ 1)Q 
2tt 



(1.16) 
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Fig. 1.5. On the left a conformal spacetime diagram showing a traveling wave in the baryon- 
photon fluid. On the right, the stationary CDM wave and the world lines of matter falling 
into the potential wells. For this wavenumber the density contrast in the baryon-photon fluid 
has undergone one-half cycle of its oscillation and is thus in phase with the Sachs-Wolfe 
effect from the CDM. This condition defines the first acoustic peak. 



the potential, so the large-scale structure (LSS) forms in the potential wells defined by the 
CDM. 

In Seljak's simple two-fluid model, there are five variables to follow: the density contrast 
in the CDM and baryons, 6 C and 5b, the velocities of these fluids v c and Vb, and the potential 
(p. The photon density contrast is (4/3)^. In Figure [T~6l the density contrasts are plotted vs. 
the scale factor for several values of the wavenumber. To make this plot the density contrasts 
were adjusted for the effect of the potential on the time, with 




Fig. 1.6. Density contrasts in the CDM and the photons for wavenumbers k = 5, 20, and 80 
(see Fig. 11.71 as a function of the scale factor relative to the scale factor and when matter 
and radiation densities were equal. The photon density contrast starts out slightly larger than 
the CDM density contrast but oscillates. 
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Fig. 1.7. Density contrasts at recombination as a function of waven umbe r k. The arrows 
on the x-axis indicate the values of k for S vs. a, as plotted in Figure fl~6l The solid curve 
shows the potential (the initial potential is always (j> = 1)> the long dashed curve curve shows 
the combined potential plus density effect on the CMB temperature, while the short dashed 
curve shows the velocity of the baryon-photon fluid. 



A e = S C + 3H J (cj)/c 2 )dt (1.17) 

and 

A 7 =<5 7 + 4// J((j)/c 2 )dt. (1.18) 

Remembering that <f> is negative when S is positive, the two terms on the right-hand side of 
the above equations cancel almost entirely at early times, leaving a small residual growing 
like a 2 prior to a eq , the scale factor when the matter density and the radiation density were 
equal. Thus these adjusted density contrasts evolve like - 1 in homogeneous universes. 

Figure ITTl shows the potential that survives to recombination and produces LSS, the po- 
tential plus density effect on the CMB temperature, and the velocity of the baryons as func- 
tion of wavenumber. Close scrutiny of the potential curve in the plot shows the baryonic 
wiggles in the LSS that may be detectable in the large redshift surveys by the 2dF and SDSS 
groups. 

A careful examination of the angular power spectrum allows several cosmological param- 
eters to be derived. The baryon to photon ratio and the dark matter to baryon density ratio 
can both be derived from the amplitudes of the first two acoustic peaks. Since the photon 
density is known precisely, the peak amplitudes determine the baryon density Ub = ftbh 2 and 
the cold dark matter density u> c = VlcDMh 2 . The matter density is given by f2,„ = Q,i, + SIcdm- 
The amplitude (Q 2 ) and spectral index n of the primordial density perturbations are also 
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Fig. 1.8. The angular size distance vs. f2 m for Ho = 60 km s" 1 Mpc" 1 and three different 

-1 12 

values of fl tot (0.9, 1, and 1.1, from top to bottom). The dashed curve shows O m ' . 



easily observed. Finally the angular scale of the peaks depends on the ratio of the angular 
size distance at recombination to the distance sound can travel before recombination. Since 
the speed of sound is close to c/y3, this sound travel distance is primarily affected by the 
age of the Universe at z = 1100. The age of the Universe goes like t oc p~ l l 2 oc 
The angular size distance is proportional to h~ l as well, so the Hubble constant cancels out. 

— 1 12 

The angular size distance is almost proportional to VL m 1 , but this relation is not quite exact. 
Figure fT~8l compares the angular size distance to flj^ 2 . One sees that a peak position that 
corresponds to il lo , = 0.95 if f2„, = 0.2 can also be fit by f2,„, = 1.1 if fi m = 1. Thus, to first 
order the peak position is a good measure of fi for . 

The CMBFAST code by Seljak & Zaldarriaga (1996) provides the ability to quickly com- 
pute the angular power spectrum Q. Typically CMBFAST runs in about 1 minute for a 
given set of cosmological parameters. However, two different groups have developed even 
faster methods to evaluate Q, Kaplinghat, Knox, & Skordis (2002) have published the Davis 
Anisotropy Shortcut (DASh), with code available for download. This program interpolates 
among precomputed Q's. Kosowsky, Milosavljevic, & Jiminez (2002) discuss combina- 
tions of the parameters that produce simple changes in the power spectrum, and also allow 
accurate and fast interpolation between Q's. These shortcuts allow the computation of a Q 
from model parameters in about 1 second. This allows the rapid computation of the likeli- 
hood of a given data set D for a set of model parameters M, L(D\M). When computing the 
likelihood for high signal-to-noise ratio observations of a small area of the sky, biases due to 
the non-Gaussian shape of the likelihood are common. This can be avoided using the offset 
log-normal form for the likelihood L(Q) advocated by Bond, Jaffe, & Knox (2000). 
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Table 1.1. Beam Size and Calibration Corrections 



Experiment 



e B (') lOO(A0)/0 lOO(AdT)/dT 



COBE 
ARCHEOPS 
BOOMERanG 
MAXIMA 



420.0 
15.0 
12.9 
10.0 
5.0 
3.0 
1.5 



-0.3 
-0.2 
10.5 
-0.9 
-0.4 
-0.5 
-1.1 



2.7 
-7.6 
0.2 
0.7 
-1.7 
-0.2 



DASI 
VSA 
CBI 



The likelihood is a probability distribution over the data, so J L(D\M)dD = 1 for any M. 
It is not a probability distribution over the models, so one should never attempt to evaluate 
J LdM. For example, one could consider the likelihood as a function of the model param- 
eters Hq in km s -1 Mpc -1 and Q, m for flat ACDM models, or one could use the parameters 
fo in seconds and D, m . For any (Ho,D, m ) there is a corresponding (/o,fi m ) that makes exactly 
the same predictions, and therefore gives the same likelihood. But the integral of the likeli- 
hood over dtodil m will be much larger than the integral of the likelihood over dHodfl,,, just 
because of the Jacobian of the transformation between the different parameter sets. 

Wright (1994) gave the example of determining the primordial power spectrum power- 
law index n, P(k) = A(k/k^) n . Marginalizing over the amplitude by integrating the likelihood 
over A gives very different results for different values of k(,. Thus, it is very unfortunate that 
Hu & Dodelson (2002) still accept integration over the likelihood. 

Instead of integrating over the likelihood one needs to define the a posteriori probability 
of the models p/(M) based on an a priori distribution pi(M) and Bayes' theorem: 



It is allowable to integrate p / over the space of models because the prior will transform when 
changing variables so as to keep the integral invariant. 

In the modeling reported here, the a priori distribution is chosen to be uniform in uj/,, uj c , 
n, Vly, £l t ot, and z„. In doing the fits, the model Q's are adjusted by a factor of exp[a + 
b£(£+ 1)] before comparison with the data. Here a is a calibration adjustment, and b is a 
beam size correction that assumes a Gaussian beam. For COBE, a is the overall amplitude 
scaling parameter instead of a calibration correction. Marginalization over the calibration 
and beam size corrections for each experiment, and the overall spectral amplitude, is done 
by maximizing the likelihood, not by integrating the likelihood. Table ll.ll gives these beam 
and calibration corrections for each experiments. All of these corrections are less than the 
quoted uncertainties for these experiments. BOOMERanG stands out in the table for having 
honestly reported its uncertainties: ±11% for the beam size and ±10% for the gain. The 
likelihood is given by 



p f (M)<xpi(M)L(D\M). 



(1-19) 
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Table 1.2. Cosmic Parameters from pre-WMAP CMB Data only 



Parameter 


Mean 


a 


Units 




0.0206 


0.0020 




UJc 


0.1343 


0.0221 




n v 


0.3947 


0.2083 




n* 


-0.0506 


0.0560 




Zri 


7.58 


3.97 




n 


0.9409 


0.0379 




Ho 


51.78 


12.26 


km s" 1 Mpc" 1 


to 


15.34 


1.60 


Gyr 


r 


0.2600 


0.0498 





+ ^/([Z°-Z?.]MZy])}, (1.20) 

i 

where j indexes over experiments, i indexes over points within each experiment, Z = ln(Q + 
Ne) in the offset log normal approach of Bond et al. (2000), and Ne is the noise bias. Since for 
COBE a is the overall normalization, a{a) is set to infinity for this term to eliminate it from 
the likelihood. The function f(x) is x 2 for small \x\ but switches to 4(|jc| - 1) when \x\ > 2. 
This downweighs outliers in the data. Most of the experiments have double tabulated their 
data. I have used both the even and odd points in my fits, but I have multiplied the er's by 
V2 to compensate. Thus, I expect to get x 2 P er degree of freedom close to 0.5 but should 
have the correct sensitivity to cosmic parameters. 

The scientific results such as the mean values and the covariance matrix of the param- 
eters can be determined by integrations over parameter space weighted by p f. Table 11.21 
shows the mean and standard deviation of the parameters determined by integrating over 
the a posteriori probability distribution of the models. The evaluation of integrals over 
multi-dimensional spaces can require a large number of function evaluations when the di- 
mensionality of the model space is large, so a Monte Carlo approach can be used. To achieve 
an accuracy of 0(e) in a Monte Carlo integration requires 0(e~ 2 ) function evaluations, while 
achieving the same accuracy with a gridding approach requires 0(e~"/ 2 ) evaluations when 
second-order methods are applied on each axis. The Monte Carlo approach is more efficient 
for more than four dimensions. When the CMB data get better, the likelihood gets more 
and more sharply peaked as a function of the parameters, so a Gaussian approximation to 
L(M) becomes more accurate, and concerns about banana-shaped confidence intervals and 
long tails in the likelihood are reduced. The Monte Carlo Markov Chain (MCMC) approach 
using the Metropolis-Hastings algorithm to generate models drawn from pj is a relatively 
fast way to evaluate these integrals (Lewis & Bridle 2002). In the MCMC, a "trial" set of 
parameters is sampled from the proposal density p,(P';P), where P is the current location 
in parameter space, and P' is the new location. Then the trial location is accepted with a 
probability given by 

P^PJP^ 
PjiP) Pt (P';P) 
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Fig. 1.9. Left: Clouds of models drawn from the a posteriori distribution based on the 
CMB data set as of 19 November 2002. The gray band shows the Big Bang nucleosynthesis 
determination of uit, (±2er) from Buries, Nollett, & Turner (2001). Right: the same set of 
models in the fi m , Ha plane. 

When a trial is accepted the Markov chain one sets P = P'. This algorithm guarantees that 
the accepted points in parameter space are sampled from the a posteriori probability distri- 
bution. 

The most common choice for the proposal density is one that depends only on the parame- 
ter change P 1 - P. If the proposal density is a symmetric function then the ratio p,(P;P')/ p,P';P) 
= 1 and A is then just the ratio of a posteriori probabilities. But the most efficient choice for 
the proposal density is p f(P) which is not a function of the parameter change, because this 
choice makes A = 1 and all trials are accepted. However, if one knew how to sample models 
from pf, why waste time calculating the likelihoods? 

Just plotting the cloud of points from MCMC gives a useful indication of the allowable 
parameter ranges that are consistent with the data. I have done some MCMC calculations 
using the DASh (Kaplinghat et al. 2002) to find the Q's. I found DASh to be user un- 
friendly and too likely to terminate instead of reporting an error for out-of-bounds parame- 
ter sets, but it was fast. Figure [l~9l shows the range of baryon and CDM densities consistent 
with the CMB data set from COBE (Bennett et al. 1996), ARCHEOPS (Amblard 2003), 
BOOMERanG (Netterfield et al. 2002), MAXIMA (Lee et al. 2001), DASI (Halverson et 
al. 2002), VSA (Scott et al. 2003), and CBI (Pearson et al. 2003), and the range of matter 
and vacuum densities consistent with these data. The Hubble constant is strongly correlated 
with position on this diagram. Figure [TTOI shows the distribution of to for models consis- 
tent with this pie-WMAP CMB data set. The relative uncertainty in to is much smaller than 
the relative uncertainty in Ho because the low-Ho models have low vacuum energy density 
(f2y), and thus low values of the product Hoto- The CMB data are giving a reasonable value 
for to without using information on the distances or ages of objects, which is an interesting 
confirmation of the Big Bang model. 

Peacock & Dodds (1994) define a shape parameter for the observed LSS power spectrum, 
F = J7,„/!exp(-2i7i). There are other slightly different definitions of V in use, but this will be 
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used consistently here. Peacock & Dodds determine T = 0.255 ± 0.017 + 0.32(n -1 - 1). The 
CMB data specify n, so the slope correction in the last term is only 0.020 ±0.013. Hence, 
the LSS power spectrum wants T = 0.275 ± 0.02. The models based only on the pre-WMAP 
CMB data give the distribution in T shown in Figure [T~TTI which is clearly consistent with 
the LSS data. 

Two examples of flat (f2, , = 1) models with equal power on all scales (n = 1), plotted on 
the pie-WMAP data set, are shown in Figure fl~T2l Both these models are acceptable fits, but 
the ACDM model is somewhat favored based on the positions of the peaks. The rise in Q at 
low I for the ACDM model is caused by the late integrated Sachs-Wolfe effect, which is due 
to the changing potential that occurs for z < 1 in this model. The potential changes because 
the density contrast stops growing when A dominates while the Universe continues to expand 
at an accelerating rate. The potential change during a photon's passage through a structure 
produces a temperature change given by AT /T = 2A(f>/c 2 (Fig. 11.131 . The factor of 2 is the 
same factor of 2 that enters into the gravitational deflection of starlight by the Sun. The effect 
should be correlated with LSS that we can see at z « 0.6. Boughn & Crittenden (2003) have 
looked for this correlation using COBE maps compared to radio source count maps from the 
NVSS, and Boughn, Crittenden, & Koehrsen (2002) have looked at the correlation of COBE 
and the X-ray background. As of now the correlation has not been seen, which is an area of 
concern for ACDM, since the (non)correlation implies = 0±0.33 with roughly Gaussian 
errors. This correlation should arise primarily from redshifts near z = 0.6, as shown in Figure 




Age [Gyr] 

Fig. 1.10. Distribution of the age of the Universe based only on the pre- WMAP CMB data. 
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Fig. 1.11. Distribution of the LSS power spectrum shape parameter V = Q m hexp(-2Q,b) 
from the pre- WMAP CMB data. 



[H4l The coming availability of LSS maps based on deep all-sky infrared surveys (Mailer 
2003) should allow a better search for this correlation. 

In addition to the late integrated Sachs-Wolfe effect from A, reionization should also 
enhance Q at low £, as would an admixture of tensor waves. Since A, r„ and T/S all 
increase Q at low £, and this increase is not seen, one has an upper limit on a weighted sum 
of all these parameters. If A is finally detected by the correlation between improved CMB 
and LSS maps, or if a substantial r„, such as the r = 0. 1 predicted by Cen (2003), is detected 
by the correlation between the Zi-mode polarization and the anisotropy (Zaldarriaga 2003), 
then one gets a greatly strengthened limit on tensor waves. 

1.3 Discussion 

The observed anisotropy of the CMB has an angular power spectrum that is in 
excellent agreement with the predictions of the ACDM model. But the CMB angular power 
spectrum is also consistent with an Einstein-de Sitter model having il m = 1 and a low value 
of Hq w 40 km s" 1 Mpc -1 . The observed lack of the expected correlation between the CMB 
and LSS due to the late integrated Sachs-Wolfe effect in ACDM slightly favors the fl m = 1 
"super Sandage" CDM model (sSCDM), which, like ACDM, is also consistent with the 
shape of the matter power spectrum P(k) and the baryon fraction in clusters of galaxies. But 
sSCDM disagrees with the actual measurements of Hq and with the supernova data for an 
accelerating Universe. Thus, ACDM is the overall best fit, but further efforts to confirm the 
CMB-LSS correlation should be encouraged. 
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Fig. 1.12. Two flat n = 1 models. One shows ACDM with fi A = 2/3. The best fit gives 
cub = 0.022 and ui c = 0.132, implying Ho = 68 km s" 1 Mpc -1 . The other fit shows J7a = 
with uj h = 0.021 and uj c = 0.196, implying H = 47 km s" 1 Mpc" 1 . 




Fig. 1.13. Fading potentials cause large-scale anisotropy correlated with LSS due to the 
late integrated Sachs-Wolfe effect. 



1.4 Conclusions 

The pre- WMAP CMB angular power spectrum assembled from multiple experi- 
ments is very well fit by a six-parameter model. Of these six parameters, the vacuum energy 
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density il v and the redshift of reionization z„ are still poorly determined from CMB data 
alone. However, the well-determined parameters either match independent determinations 
or the expectations from inflation: 

• The baryon density u>b is determined to 10% and agrees with the value from Big Bang 

nucleosynthesis. 

• The age of the Universe is determined to 11% and agrees with determinations from white 

dwarf cooling (Rich 2003), main sequence turnoffs, and radioactive decay. 

• The predicted shape of the LSS power spectrum P(k) agrees with the observed shape. 

• The curvature is determined to 4% and agrees with the expected value from inflation. 

• The spectral index n is determined to 4% and agrees with the expected value from inflation. 

The angular power spectrum of the CMB can be computed using well-understood physics 
and linear perturbation theory. The current data set agrees with the predictions of inflation 
happening less than 1 picosecond after the Big Bang, the observations of light isotope abun- 
dances from the first three minutes after the Big Bang, and the observations of LSS in the 
current Universe. The inflationary scenario and the hot Big Bang model appear to be solidly 
based on confirmed quantitative predictions. 

The greatly improved CMB data expected from WMAP, and later Planck, should dramat- 
ically improve our knowledge of the Universe. 
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